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The impact of the oblateness of Regulus on the motion of 
its companion 
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Abstract The fast spinning B-star Regulus has re- 
cently been found to be orbited by a fainter com- 
panion in a close circular path with orbital period 
Pb = 40.11(2) d. Being its equatorial radius R c 32% 
larger than the polar one R p , Regulus possesses a re- 
markable quadrupole mass moment Q. We investigate 
the effects of Q on the orbital period Pb of its com- 
panion in order to see if they are measurable, given the 
present-day level of accuracy in measuring Py,. Con- 
versely, we will look for deviations from the third Ke- 
pler law, attributed to the quadrupole mass moment 
Q of Regulus, to constrain the ratio 7 = m/M of the 
system's masses. The impact of Q on the orbital pe- 
riod is analytically worked out with a straightforward 
perturbative approach. The resulting correction P® is 
compared to other competing dynamical effects. P^ 
and the Keplerian period P Ke P are expressed in terms 
of the phcnomcnologically determined system's param- 
eters; 7 is treated as an unknown. P® is compared to 
the observational accuracy in measuring the orbital pe- 
riod 5Pb = 0.02 d and to the systematic uncertainty 
<5(P Kep ) due to the errors in the system's parameters 
entering it. The discrepancy AP = |Pb — P Kc p| is ex- 
amined in order to sec for which values of 7 it becomes 
statistically significant. The physical meaning of the 
obtained range of values for 7 is discussed in terms of 
Q. P® is larger than <5Pb but still smaller than the 
systematic uncertainty in P k °p by two orders of mag- 
nitude. The major sources of bias are the velocity semi- 
amplitude K of the motion of the primary and its mass 
M, Assuming edge-on configuration, i.e. i = 90 deg, if 
7 > 0.096 Q would be positive, i.e. Regulus would be 
prolate, contrary to the observations. If 7 < 0.078 Q 
would be negative, but its magnitude would be one-two 

Lorenzo Iorio 

INFN-Sezione di Pisa. Permanent address for correspondence: 
Viale Unita di Italia 68, 70125, Bari (BA), Italy. 



orders of magnitude larger than the approximate esti- 
mate Q « M(Rl-Rl) = -2.4±0.5x 10 49 kg m 2 . Reg- 
ulus is the first extrasolar binary system in which the or- 
bital effects of the asphericity of the primary are larger 
than the observational sensitivity; moreover, no other 
competing aliasing orbital effects are present. Thus, it 
is desirable that it will become the object of future in- 
tensive observational campaigns in order to reduce the 
systematic uncertainty due to the system's parameters 
below the mcasurability threshold. 

Keywords Stars: individual: Regulus; Stars: bina- 
ries: close; Stars: fundamental parameters Gravitation; 
Celestial mechanics 



1 Introduction 

Regulus (a Leo, H R 3982, HP 8790 1) is a nearby 
(d = 24.3 ± 0.2 pc (Ivan Leeuwen 2007)) intermediate 
mass star of spec tral class B (jjohnson fc Morganlll953 : 
iGrav et al.ll2003h and has the peculiarity of being ani- 
mated by a very fast rotation. Indeed, combined inter 



ferometric and spectroscopic studies (jMcAlister et al 



20051 ) have shown that its equatorial radius R B is 32% 
larger than the polar one P p , and its rotation period 
is 15.9 hr, which corresponds to an equatorial rotation 
speed that is 86% of the critical breakup velocity. Other 
fast spinning stars whose oblateness has been interfer- 
ometrically measured are Altair (a Aquilae, HP 187642 

( 

( 



van Bell e et alJl200lft1 and Ache rnar (a Eri, HP 10144 



Pomiciano de Souza et al.ll2003l )) 



In principle, a companion orbiting not too far from 
such highly deformed stars would experience rele- 
vant orbital effects induced by the primary's oblate- 
ness which could be used to dynamically put con- 
straints cither on such an important stellar physi- 
cal parameter or on the orbital/physical parameters 



Table 1 Relevant parameters of the single-line spec- 
troscopic bin ary system of Regu lus. M is the pri- 
marv' s mass JMcAlister et al.l I2005T ). Ph is the orbital pe- 
riod (IGies et al.l 120081 ). K is the velocity semiamplitudc 
(JGies et al.l 120081), R e and R v are the equatorial and po- 
lar radii ( McAlister et al]|2005f): the ec centricity e has been 
assumed equal to zero ( Gies et al.ll2008l ). Numbers in paren- 
theses give the error in the last digit quoted. 

M (M ) P h (d) K [*f) R c (Rq) R V ~J^T 
3.4(2) 40.11(2) 7.7(3) 4.16(8) 3.14(6) 



of the system. In fact, the Be star Achernar re- 
cently turned out to be orbited by a f aint companion 
1 Kervella fc Domiciano de Souza 2007 ) w hich should 
be a A1V-A3V star (JKervella et alj l2008h . but since 
its orbital period should amount to about 15 yr it has 
not yet been possible to obtain its orbital parameters. 
No companions are known for Altair. The situation is 
more favorable for Regulus since spectroscopic obser- 
vations over the last few years have demonstrated that 
it has a clo sqj companion i n a circular orbit described 



in 40.11 d dGies et al.ll20ol) 



In this paper we will investigate the possibility of de- 
tecting some dynamical orbital effects induced by the 
oblatcncss of Regulus on its companion and will use 
them to put some constraints on the mass of the sec- 
ondary. 



2 The impact of Q on the orbital period and 
its measurability 

The relevant physical and orbital parameters of the 
Regulus binary system are in Table [TJ Concerning the 



McAlister et al 



inclina tion i to the plane of the sky, _ 
( 2005T ) showed that the best results for their fits are ob 
tained by choosing i = 90 deg; in the following we will 
keep it fixed to different values close to the edge-on con- 
figuration. The estimates of stellar radii and mass are 



almo st independent of the choice of i (jMcAlister et al 
2005h . 

Given that the observable quantity at our disposal 
is the orbital period, we will start by investigating the 
impact of the primary's oblateness on it. By assuming 
axial symmetry about the z axis and reflection symme- 
try about the equator assumed as reference {xy} plane, 
the external gravitational potential U can be written as 



1 In fact, Regulus has a known wide companio n a Leo B at a 
separ ation of « 175" which is itself a binary JMcAlister et al.l 
[2005J) , but it has a far too great separation to have ever interacted 
directly with Regulus. 



( Laarakkers fc Poissonlll999[ ) 

U = U Q + U Q , 
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GM 
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(1) 



(2) 



in which 9 is the eolatitude angle and Q = Q + q ss Q; 
indeed, the quadrupole mass moment is proportional 
to the square of the stellar angular rotation frequency, 
so that we will neglect the oblateness q, if any, of the 
secondary with respect to that of Regulus. The rela- 
tive acceleration A = — VC/ due to the gravitational 
potential of eq. (J2J) is, in spherical coordinates 



n-r — g r , 



A, 
A, 
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which yields 
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(4) 



A v = 0. 



We will now make the simplifying assumption that 
the orbital ang ular momentum and the spin of Regu - 
lus are aligned (JMcAlister et al.ll2005l : lGies et al.ll2008L 



Thus, Ag = A v = and only the equation for the radial 
acceleration survives in eq. (U]) as 



(5) 




(6) 



Following straightfor wardly a standard perturbative 
approach ( Ioriol I2008T ) is it possible to use eq. 



obtain 

p — pKcp , pQ 

with 
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Kep _ 
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to 



(7) 



(8) 
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Fig. 1 

Keplerian period for M 



Quadrupole mass moment correction P y to the 
3.4M , Q » M(R 2 p -R 2 C ) = 
(—2.4 ± 0.5) x 10 49 kg m 2 and different values of the incli- 
nation i close to the edge-on configuration. It is larger than 
the accuracy in the determined orbital period <SPb = 0.02 d. 



Fig. 2 Total uncertainty 8(P Kep ) in the Keplerian period 
due to the errors in the system's parameters for M = 3.4Mq 
and different values of the inclination i close to the edge-on 
configuration. It is larger than P® by about two orders of 
magnitude. 
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M = M + m = M(l + 7), 7 
and 
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(10) 



(11) 



where a is the relative semimajor axis and om is the 
barycentric semimajor axis of the motion of the pri- 
mary. 



Since (JLaarakkers fc Poissonlll999l ) 



Q = Ix-h~ M{Rl - Rl) = -2.4 ± 0.5 x 10 49 kg m 2 , 

(12) 

where I x and I z are the equatorial and polar moments 
of inertia, respectively, eq. © allows us to investigate 
if the contribution of the primary's oblateness to the or- 
bital period is large enough to be detected, in principle, 
given the present-day level of accuracy in determining 
Pb- Figure[T]and Tablc[T]show that this is just the case. 
Here and in the following we will treat 7 as an unknown 
quantity. This is an important result because it opens, 
in principle, the possibility of inferring or constraining 
Q from the dynamics of the binary system. Moreover, 
it is the first time that such an opportunity occurs for 
an extrasolar binary system; indeed, we will show that 



other competing dynamical effects, which in other sys- 
tems may overcome those due to Q, in this case are 
negligible. 

To do that it would be necessary compute the Ke- 
plerian orbital period P Ko p and subtract it from the 
measured one Pb; to this aim, let us note that the pa- 
rameters entering eq. ((SJ have been determined inde- 
pendently of the third Kepler law itself, so that it does 
make sense to compare eq. © to the phenomenologi- 
cally determined orbital period Pb- However, it is not 
sufficient that 5P h < P Q ; is the uncertainty <5(P Kcp ) 
due to the errors in the system's parameters smaller 
than P^ as well? Unfortunately, the answer is still 
negative by about two orders of magnitude, as shown 
by Figure [U 

It may be interesting to see what is the impact of 
each system's parameter in determining <5(P Kcp ) in or- 
der to have some hints about the future possibilities 
and also to drive researches towards the improvement 
of the parameters which turn out to be the most ef- 
fective in corrupting the measurement of P®. Figure 
[3] shows that the most relevant sources of errors are K 
and M. For details concerning ho w K and M and their 
errors have been deter mined see (IGies et alj|2008l) and 
references therein, and (JMcAlister et alj|200a ). respec- 
tively. 

In principle, there is also a general relativistic con- 
tribution to the or bital period which should be taken 
into account. It is ( Damour fc Deruellelll986l ) 



F*«£VGM5=l(l+7)< 
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Fig. 3 Uncertainties in the Keplerian period due to the 
errors in the system's parameters M, K, p, for M = 3.4Mq 
and i — 90 deg. 



However, it can be safely neglected because it is smaller 
than SPh by four orders of magnitude, being of the order 
of 1CT 6 d. 

Finally, let us note that, in view of the large sepa- 
ration between the two stars no tidal effects come into 
play. This can be seen by comparing the tidal acceler- 
ation 



A mal = ^^ = 2GM 7 4 i? c 



2-7T sin i 



(l+ 7 )KP b 



to the centrifugal acceleration 



(14) 



(15) 



where ft is the rotational frequency of Regulus; it turns 
out that j4tidai/^4con "C 1 being of the order of 10~ 9 . 



3 Constraints on the secondary's mass from 
deviations from the third Kepler law 
attributed to Q 

In this Section we will show how useful constraints on 
the mass m of the Regulus' companion can be inferred 
by looking for deviations from the third Kepler law in 
terms of Q; indeed, in the previous Section we showed 
that no other effects are relevant in the orbital dynamics 
of such a binary system. 

Let us start by evaluating the discrepancy AP = 
\P\, — P Kcp \ between the measured orbital period Pb 
and the computed Keplerian one P Ke P as a function 
of 7 for different values of the inclination close to the 
edge-on configuration: Figure [3] shows that, for i = 90 
deg, AP is significant at 1 — a level outside the range 
0.081 < 7 < 0.093; for i = 75 deg it occurs outside 



Fig. 4 Ratio AP/5(AP) of the discrepancy AP = |P b - 
P Kep | between the measured orbital period Pb and the Kep- 
lerian one P Kcp to its uncertainty as a function of 7 = m/M 
for M — 3.4M© and different values of the inclination i close 
to the edge-on configuration. It is significant at 1 — a level 
outside the range 0.081 < 7 < 0.093 (i = 90 deg). 



0.085 < 7 < 0.096. Note that such bounds are conser- 
vative because the uncertainty in AP has been evalu- 
ated as S(AP) < SP h + S(P Kc p). For the sake of sim- 
plicity, let us consider the case i = 90 deg and see if 
statistically significant deviations from the third Kepler 
law for 7 > 0.09, interpreted as due to Q, are physically 
meaningful. Figure \E\ clearly shows that the answer is 
negative: indeed, for 7 > 0.096 the quadrupole mass 
moment Q of the primary would become positive, i.e. 
the star would be prolate, which is contradicted by the 
observations. Thus, we can conclude that the faint com- 
panion of Regulus cannot have a mass larger than about 
0.30 Mm , contrary to what can be found in (JGies et al 
20081) . 



An inspection of Figure [5] tell us that for 7 < 0.078, 
i.e. m < 0.26M Q , Q would be a definite negative quan- 
tity, in according with the observations. However, it 
must be noted that, in this case, \Q\ would be one- 
two orders of magnitude larger than the estimate of eq. 
(fl"2"l) . It would be difficult to consider such a large dis- 
crepancy as acceptable. Thus, we provisionally consider 
7 < 0.078 not likely. 

The simplest solution consists in assuming for 7 the 
range which still makes AP compatible with zero con- 
sidering also that the value of eq. (fT2")l for Q would not 
be excluded by such values of 7, as shown by Figure [5l 



4 Conclusions 

In this paper we investigated the impact of the huge 
quadrupole mass moment Q of Regulus on the orbital 
period of its faint companion. 
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Fig. 5 Allowed values for the quadrupole mass moment 
Q as a function of 7 = m/M for M = 3.4M Q and i = 90 
deg. For 7 > 0.096 Q becomes positive, i.e. the star would 
be prolate, contrary to the observational evidence. For 7 < 
0.078 Q is negative, but one-two orders of magnitude larger 
than M(B% --Re). 



It turns out that, for a reasonable estimate of its 
value, Q induces a correction to the Keplerian period 
which could, in principle, be measured because it is 
larger than the error in phcnomcnologically measur- 
ing the orbital period Pb- However, according to the 
present-day level of knowledge of the system's param- 
eters, the uncertainty in the Keplerian period is larger 
than the correction due to Q by about two orders of 
magnitude; the most important sources of errors are 
the velocity semiamplitude K and the mass M of Reg- 
ulus. 

An analysis of deviations from the third Kepler law 
as due to Q for the edge-on orbital configuration and 
M = 3.4(2)M Q has shown that values of the ratio 7 of 
the secondary's mass to the primary's one larger than 
about 0.096 are to be ruled out because they would 
yield a positive quadrupole mass moment for Regulus, 
i.e. it would be prolate, contrary to the observations. 
On the other hand, for 7 < 0.078 Q would be negative, 
but with an unlikely large value. 

Further observational campaigns with different tech- 
niques would be of great significance in the attempt of 
reducing the overall uncertainty in the system's param- 
eters down to the level required for detecting the effect 
of the shape of Regulus. 
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